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ABSTRACT 

We present a detailed study of the evolution of GMCs in a galactic disc simulation. We follow 
individual GMCs (defined in our simulations by a total column density criterion), including 
their level of star formation, from their formation to dispersal. We find the evolution of GMCs 
is highly complex. GMCs often form from a combination of smaller clouds and ambient ISM, 
and similarly disperse by splitting into a number of smaller clouds and ambient ISM. However 
some clouds emerge as the result of the disruption of a more massive GMC, rather than from 
the assembly of smaller clouds. Likewise in some cases, clouds accrete onto more massive 
clouds rather than disperse. Because of the difficulty of determining a precursor or successor 
of a given GMC, determining GMC histories and lifetimes is highly non-trivial. Using a 
definition relating to the continuous evolution of a cloud, we obtain lifetimes typically of 4- 
25 Myr for > 10^ M© GMCs, over which time the star formation efficiency is about 1 %. 
We also relate the lifetime of GMCs to their crossing time. We find that the crossing time 
is a reasonable measure of the actual lifetime of the cloud, although there is considerable 
scatter. The scatter is found to be unavoidable because of the complex and varied shapes and 
dynamics of the clouds. We study cloud dispersal in detail and find both stellar feedback and 
shear contribute to cloud disruption. We also demonstrate that GMCs do not behave as ridge 
clouds, rather massive spiral arm GMCs evolve into smaller clouds in inter-arm spurs. 
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1 INTRODUCTION 

Observational studies of molecular clouds are limited by providing 
only an instantaneous picture of where clouds lie and what their 
properties are. As such we have a restricted understanding of the 
evolution of molecular clouds, for example of how they form, what 
they form from, and for how long they undergo star formation. Nu- 
merical simulations offer one way to study cloud evolution. How- 
ever even with simulations, studying the evolution of clouds is still 
challenging, with problems such as how to define a clump or cloud 
(see e.g. |Pineda et al.|20 09; Reid et al.|2010l [Hughes et al.,|2013| ), 
and because the constituent gas of clouds may change on relatively 
short timescales (e.g. Myrs, Dobbs et al.|2011b) ). 

Observations of Galactic GMCs show that some, such as 
Orion, are actively forming stars, whereas others, such as Califor- 
nia and Maddalena's cloud are much more quiescent ( Lada et al. 
|2009[pegeath et al. 20 09). Nevertheless all nearby GMCs appear 
to contain at least some star formation, seemingly since all con- 
tain some fraction of very dense (^ 10^ cm~^) gas ( Froebrich 



|& Rowles 2010, Lada et al.|,2010 K One implication of these ob 
servations, and the generally small age spread of stars in GMCs 
pfremov & Elmegreen|[T998l [Mayne & Naylor|[2008l [Sung &] 
|Bessell||2010| [Jeffries et al.||201li ), is that the timescale for star 



formation is very short ( [Larson|1981| [Elmegreen|2000| [Hartmann| 
jet al.|200"T] ). Likewise this implies a short formation time for GMCs 
themselves ( [Pringle et al.|200l| ). To account for these short forma- 
tion times, there have been two main possibilities - either diffuse 
gas is collected together in a colliding flow (Audit & Hennebellej 



2005| jHeitsch et al. [[20061 jVazquez-Semadeni al.|[2007[ [20101 



Inoue & Inutsuka[[20 12^ on a Myr timescale, or gas is relatively 
dense prior to becoming a GMC ( Pringle et al. 200T]). In [Dobbs] 
jet al.H 2012 ) we considered the gas which goes to form the global 
population of GMCs at a given timeframe, and found that gas is 
typically overly dense for tens of Myr prior to GMC formation, 
supporting the latter scenario. Observing the precursors of GMCs 
is difficult and it is not known whether they are dense clouds of HI, 
diffuse HI, or other molecular clouds. In M33 and the LMC, there 
have been observations of GMCs that do not contain massive star 
formation ( [Kawamura et al.[[2009] jGratier et aTl|2012| ), but there 
may well be unobserved low mass star formation in these clouds, 
similar to the California cloud. 

Numerical simulations have also investigated the physical 
mechanisms for GMC formation (Dobbs et al. 2006 , Dobbs 2008', 
Wada 2008 , IKim & Ostriker 2002, Kim et al. 2003 , Tasker & Tan] 



2009}. Inp)obbs (2008) we showed that at low surface densities 
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(< 10 M© pc^) GMCs form from the collision or agglomeration 
of smaller clouds. As the surface density increases, gravitational 
attraction between the clouds enhances collisions ( [Kwan & Valdesj 
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|1987| ) and gravitational instabilities in the gas play a more impor- 
tant role ( |Kim & Ostriker|2002||2006] ). This picture of agglomer- 
ation and self gravity is largely supported by recent observations 
of the GMCs and the properties of GMCs in M51 ( Colombo et aU] 
|20T3l ). 

The lifetimes of GMCs are not well known either observa- 
tionally or theoretically (although as we point out in this paper the 
definition of a GMC lifetime is not clear either), although many 
observations point towards a lifetime of 20 — 30 Myr. The first 
observational estimates based on the presence of CO surrounding 
clusters of different ages, suggested that the ages of Galactic GMCs 
are - 20 Myr ( |Bash et al. 1977 , Leisawitz et al. 1989 , Blitz 1993 ). 
More recent results have estimated cloud lifetimes in the LMC and 
M33 to be 20-30 Myr (Kawamura et al.|2009[piura et aLpOTI] ). 
[Murray! ( |2011| ) deduced similar lifetimes based on the estimated 
time for feedback to disrupt observed GMCs. There have however 
been suggestions of longer lifetimes, for example to explain the 
presence of inter-arm GMCs in the Milky Way ( Scovill e etal.|1979| ) 
and other galaxies such as M51 ( ^Koda et al. 2009j . Over a cloud's 
lifetime some fraction of the cloud's mass it turned into stars, rep- 
resenting the star formation efficiency. This fraction is typically a 
few per cent and is thought to range from < 1% to 20% ( [Evans | 
[et al. [2009] [Murray [2011] ). The relatively low value is thought to be 
due to stellar feedback, which changes the dynamics and disrupts 
GMCs ( [Dobbs et al.|2011a|rV azquez-Semadeni et al. 2010 ) and/or 
magnetic fields (P rice & Bate|2008nPadoan & Nordlund|2011| . 

The primary method of disruption of GMCs is also not well 
known. GMCs are subject to external processes, such as shear and 
turbulence, which can disperse the cloud, as well as internal pro- 
cesses of star formation and stellar feedback. Recently ^ Dobbs et aT] 
( [201 lb| ) argued that many or most GMCs in galaxies are unbound. 
For these, internal processes such as stellar feedback are not nec- 
essary to disrupt the cloud. For bound GMCs, and for bound sub- 
structures within clouds, disruption requires energy input from stel- 
lar feedback processes such as jets, winds, HII regions and super- 
novae. 

In a previous paper ( [Dobbs etal.|2012[ ) we investigated the gas 
density of material before and after it forms molecular clouds. We 
found that gas is likely to be overly dense compared to the aver- 
age densities of the ISM, at least 30 Myr prior to GMC formation. 
In particular, gas which reaches moderate densities of ^ 10 cm~^ 
persists for tens of Myr, even though dense 100 cm~^) gas is 
only very short-lived. Furthermore we found some examples where 
clouds were never completely dispersed into the ambient medium. 
Here we extend these ideas, but to show examples of individual 
clouds in detail. We follow their formation, dispersal and star for- 
mation history. We focus on a representative 10^ M© cloud, but 
also describe the evolution of other 10^ — 10^ M© clouds, and 
show some overall statistics for > 10^ M© clouds. We also vary 
our criteria for defining a cloud to check the reliability of our re- 
sults. 



2 CALCULATIONS 

We performed a calculation with a 2 armed spiral potential, heat- 
ing and cooling, self gravity and stellar feedback. The implemen- 
tation of all these processes has been described in previous papers. 
The potential used comprises a logarithmic component ( Bi nney &[ 
[Tremaine]l987 ) and a spiral component ( i Cox & G6meZ|200 2l), and 
is the same as that given in [Dobbs et al.| ( [2006| ), except in this pa- 
per we choose N — 2 rather than N — A. We changed io N — 2 



simply because for many roughly symmetric galaxies, the number 
of arms is 2 rather than 4 (the choice of the latter previously was 
also influenced partly by the supposed structure of our Galaxy (see 
[Cox & G6mez|2002 and Vallee 200 5 )). 

The cooling and heating of the ISM we use is taken from 
Glover & Mac Low] (12007^ and described more fully in D obbs] 
et al.[ ( [2008| ). Cooling is switched off if the temperature of a parti- 
cle drops below 50 K. Also as in these previous papers we include 
time-dependent H2 formation, described briefly as follows. Each 
particle is assigned a fraction of H2 (from to 1) and the rate of 
change of n{H2) is calculated according to the prescription given 
in the appendix of Bergi n et al.[ ( [2004| ), and [Dobbs et al.[ ( [2006 ]). 
This prescription includes the rate of formation of H2 on grains 
(see e.g. Hollenbach & "Salpeter| 197 1 [ [Jura| 1974 j and a destructive 
photodissociation term ( [Draine & Bertoldi[1996| ). For the photodis- 
sociation term, we used a simple approximation of estimating the 
column density from the local density multiplied by a scale-length 
of 35 pc, as described in [Dobbs et al.[ ( [2008| . 

We implement stellar feedback as described in [Dobbs et al.[ 
( [201 la| ), but also provide the main details here. Although some re- 
cent work has considered in detail different mechanisms such as 
ionisation, winds and supernovae ( Hopkins et al.|201 l[p\gertz et aT] 
[2012|[Kiin et al.|20T2| , we simply insert an amount of energy where 
star formation occurs. We insert an amount of energy according to 

eM{H2) ,^51 
160 

where e is a constant which we take to be 0.05, M(H2) is the mass of 
molecular hydrogen in a region of gas (dependent on the smooth- 
ing length at a particular resolution), the division by 160 comes 
from our choice of IMF, and we multiply by 10^^ ergs assuming 
this is the energy released from one supernovae. For our resolution, 
each feedback event typically corresponds to 1 or 2 supernovae. 
However we note that our scheme could count collectively for all 
stellar feedback. The energy associated with processes such as jets 
is much smaller, whilst that for stellar winds and radiation pressure 
are comparable to supernovae (Mac Low & Klessen 2004} [Agertz[ 
[et al.[2012] ), hence the net amount will not be greatly different from 
the 10 ergs we assume. We deposit the energy as half kinetic and 
half thermal energy. 

The simulation presented here uses 8 million particles, each 
particle having a mass of 312.5 M©. We also present a 1 million 
particle simulation in the appendix as a resolution test. The mean 
surface density of the galaxy is 8 M© pc~^. 



3 RESULTS: EVOLUTION OF GALAXY 

In Figure 1 we show the column density for our galactic disc at a 
time of 250 Myr. The galaxy exhibits two clear spiral arms in the 
central parts, though the arms become weaker at larger radii. The 
difference in velocities across the shock drop substantially beyond 
^ 6 kpc reflecting that, although the depth of the potential (the 
spiral perturbation is of order 0.1 % of the logarithmic component) 
does not alter significantly, the potential minimum is much broader 
and the gas does not undergo a strong shock. 

The gas in the disc is clumpy, with the most massive, dense 
clumps or clouds situated in or near the spiral arms. These massive 
clouds form by a process of agglomeration of smaller clumps in 
the spiral arms, aided by self gravity. Clouds in the spiral arms are 
then sheared in to inter- arm spurs, producing a structure remarkably 
similar to observed grand design galaxies such as M5 1 and M74. 
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Figure 1. The column density for the disc in our simulation is shown at a 
time of 250 Myr. The red box indicates the region shown in Figure 3 (middle 
panel) which contains Cloud380. 

From 150 Myr onwards, the gas lies in approximately one third of 
the cold (< 150 K), intermediate and warm (> 7000 K) phases of 
the ISM respectively. A small fraction of the gas, ^ 1 %, lies in the 
hot (> 10^ K) phase. We ran our calculation until a time of 300 
Myr, and we carry out most of our analysis at, or around the time 
of 250 Myr. 



4 DEFINING AND TRACING CLOUDS 

In this section, we describe how we select clouds in our simulation, 
how we trace them over time, and how we estimate the lifetime of 
a cloud. We demonstrate that the clouds do not behave as 'ridge' 
clouds, so it appropriate to follow them as a collection of parti- 
cles. We also discuss how cloud can be dispersed, both in terms of 
physical processes and with regards to the limitations of numerical 
simulations. 



the properties we obtain are similar to the properties of clouds pre- 
sented in pobbs et al.|pOlTa| , although in that calculation we only 
used 1 million particles, and adopted a 4 arm spiral potential. We 
note also that the properties of GMCs in our main, 8 million parti- 
cle simulation, extend to lower mass (10^ M©) clouds compared to 
IDobbs etal.| ( |2011a| ). 



4.3 Tracing the evolution of clouds 

Since we are using a Lagrangian code, we are able to trace the 
constituent particles in a given cloud to earlier and to later times. 
However, when a cloud is defined at a particular moment in time 
(for example by using a clump-finding algorithm) it is only really 
defined at that instant. At the next instant (or at a previous instant) 
what appears to be the same cloud will in fact be a new cloud with 
(slightly) different constituent gas and different boundaries. At later 
times, the cloud may break up into one or more clouds - some of 
the original cloud gas will be split among many other clouds, and 
some may have returned to the intercloud medium. Similarly, if we 
try to trace back the cloud to earlier times, the gas in the cloud will 
have come from a number of other clouds and/or directly from the 
intercloud medium. Thus the concept of the evolutionary history 
and lifetime of a cloud is not necessarily a particularly useful one. 

With this in mind, we set about tracing the evolution of a cloud 
in the following manner: (i) at a particular time To (here chosen to 
be 250 Myr) we consider a particular cloud A defined by the clump 
finding algorithm, and nominate the constituent particles of cloud 
A at that time as 'original cloud particles', then (ii) at other times 
t / To we use the same clump finding algorithm and list all those 
clouds which contain more than 20 (i.e. > 625OM0) of the original 
cloud particles (Recall that by definition a cloud must contain more 
than 30 particles.). Then at each time t we can identify that cloud 
which contains the greatest number of original cloud particles, and 
we nominate that cloud as the continuation of cloud A in time. 

In this way we are able to establish the evolution of a cloud in 
a Lagrangian sense, that is we define the evolving cloud in terms of 
the movement of a set of original cloud particles. 



4.1 Cloud selection 

We used a simple clump-finding algorithm which, for most of our 
analysis, selects clumps according to a surface density threshold 
criterion of 12th = 100 Mq pc~^. The galaxy is divided into a 
Cartesian grid, with cells of size 10 pc, and all adjacent cells which 
exceed the density threshold are grouped together and classified 
as a cloud (cells which are only diagonally touching are not in- 
cluded). We also require that clouds contain at least 30 particles. 
At 250 Myr, this gave 1562 clouds with masses in the range 10^ 
M0 to 2 X 10^ M0. We refer to any material which is not in a 
cloud as 'ambient ISM'. In Section 5.3, and the appendix, we also 
show results where we vary the surface density threshold (referred 
to hereafter as St/^) used in our clump-finding algorithm. 

4.2 Properties of clouds 

In Appendix A, we show the mass spectrum, rotation and virial 
parameters of the clouds we detect. In addition, we also compare 
results with different surface density criteria, and results from a 
lower resolution simulation with only 1 million particles. Overall 



4.3.1 Molecular clouds are not 'ridge clouds' 

As an aside here we note that this is not necessarily the only, 
nor the best, definition of cloud evolution, depending on circum- 
stances. For example if the molecular cloud were the astronomical 
equivalent of a 'ridge cloud' (that is, a rain cloud that forms near 
a mountain ridge as air is forced over it), perhaps being caused 
as gas is compressed by passing through a spiral arm, then the 
cloud would remain roughly in one place as its constituent parti- 
cles passed through it. 

To show that this would not be a sensible procedure in this 
case we consider Cloud380 described in more detail in the next 
Section. At time To we define the volume V which coincides with 
Cloud380. We then plot the mass of gas Mv{t) which lies within 
this volume (in the co-rotating frame) in Figure 2, and compare 
this with the mass of original cloud particles Mvo (t) within V at 
each time t. Were Cloud 380 a ridge cloud, then Mv{t) would re- 
main roughly constant, while the 'original cloud particles' Mvo{t) 
would flow through the volume. We can see from Figure 2 that My 
and Mvo track each other quite well. This suggests that Cloud 380 
is not a ridge cloud, and that tracking its evolution by following its 
constituent particles does make sense. 



© 2012 RAS, MNRAS 000,[l]{T5] 



4 Dobbs 




230 235 240 245 250 255 260 
Time (Myr) 

Figure 2. The change in mass of gas in a volume V is shown versus time. 
V includes, and is slightly larger than Cloud380 (see Section 4.4). The two 
lines show the total mass of gas in V versus time {My (t) - see text, blue 
dashed) and the mass of the original particles of V at time Tq = 250 Myr 
(MyQ (t), red solid) versus time. Since the evolution of both is very similar, 
we conclude that cloud evolution is essentially Lagrangian and the cloud is 
not a ridge cloud. 



4.3.2 Dispersal of clouds 

As we show later, the dispersal of clouds in the simulations is 
highly complex. There are several physical mechanisms by which 
clouds can be dispersed, i) the simple unboundedness of the cloud, 
ii) shear, and iii) stellar feedback, or some combination of these. 
The complex structure of many clouds means it is easier to split 
clouds apart by these processes, e.g. if a cloud is a long filament, 
it is relatively easy to split the cloud into two if there is a feed- 
back event, the cloud is subject to shear, or generally one part of 
the cloud is experiencing different large scale motions to the rest. 

There are also a couple of artificial ways that clouds can dis- 
perse (or form), due to the resolution of our simulation and method 
for finding clouds. A cloud can disappear completely simply by 
falling under the threshold used to determine clouds. For more mas- 
sive clouds though it is unlikely that the whole cloud will disappear 
simultaneously this way. Also clouds form and disappear instanta- 
neously if the number of particles changes between 29 and 30. As 
we restrict our studies to > 10^ M© clouds with > 320 particles, 
the artificial dispersal (and formation) of these clouds is likely less 
important, but worth bearing in mind. 



4.3.3 Cloud lifetimes 

Because, as we show in the next sections, it is often difficult to 
establish the precursor or successor of a given cloud, it is very dif- 
ficult to assign a lifetime to a given cloud. However we can use the 
continuation of a given cloud to assign its lifetime. We compute a 
lifetime as the time at which the continuation of a given cloud is 
at least half the mass of that cloud at To (see Sections 5.1.5 and 
5.3.1). This is a somewhat arbitrary definition, but does allow us to 
quantitatively compare clouds. 



5 RESULTS 

To start, we consider in detail one particular cloud in our simula- 
tion, Cloud380. In Section 5.2. we consider several other massive 
clouds in the simulation, then in Section 5.3 we consider statistics 
of all clouds over 10^ M0. 

5.1 Example 1: Cloud380 

We consider first one of the more massive clouds found by the 
clump finding algorithm at time To = 250 Myr and name it Cloud 
380 (being number 380 on the algorithms list). It is located in a 
spiral arm 3 kpc from the centre of the galaxy, and has a mass of 
2 X 10^ M©, corresponding to 6386 particles. In Figure 3 we show 
the appearance of the cloud at time t = To = 250 Myr (centre 
panel), and show the evolution of the cloud material over a period 
230 Myr < t < 270 Myr. The mass of the continuation of Cloud 
380 is shown in Section 5.1.5. In the vertical plane. Cloud 380 ap- 
pears as one coherent structure, centred close to the midplane of the 
galaxy (see also |Dobbs et al.|2011a| for examples of edge on plots 
of simulated galaxies). 

5.1.1 Cloud dispersal: evolution at times t > Tq 

If we start with the centre panel of Figure 3 and move forward in 
time, we find that after 20 Myr only ten per cent of the original 
cloud is still in clouds - the rest has been returned to the ambient 
ISM. In each of the later panels we delineate approximately the 
area covered by the original cloud particles. The major successor 
cloud is shown in red, but there are also other clouds (shown in 
different colours) which contain a substantial number (> 20) of 
the original cloud particles. The numbers by some of these clouds 
shows what fraction of each of those clouds is made up of original 
cloud particles. Thus we see that after To Cloud380 fragments into 
multiple smaller clouds, although there is still one massive cloud 
evident. The cloud disperses through a combination of shear and 
feedback (We discuss the role, and relative importance, of these 
processes below.). 

5.1.2 Clouds, spurs and the effects of stellar feedback 

Until To = 250 Myr, Cloud380 (or rather the progenitor of 
Cloud380) experiences net growth as it travels along the spiral arm. 
However after 250 Myr, there is not much gas in the spiral arm in 
the vicinity of Cloud380, and Cloud380 then leaves the spiral arm, 
hence its mass decreases. Not long after To, for 255 Myr < t < 
270 Myr, the whole of the successor of Cloud380 lies in a 'spur' 
feature downstream of the arm. At times from 255 to 265 Myr, the 
successor of Cloud380 is the largest cloud within a 1 kpc radius by 
a factor of a few. In particular there are no comparable mass clouds 
in the spiral arm. After 270 Myr, material begins to build up in the 
spiral arm again. At 270 Myr, the most massive GMCs in the arm 
and spur are comparable (^-^ 2 x 10^ M©). 

This picture can be compared with the results of |Schinnerer et| 
|al.,|p013} where they study clouds and spurs in a region of M51. 
They find that the most massive clouds are preferentially located on 
the downstream side of the arm. For their particular section of a spi- 
ral arm, Schinnerer et al.,| ( |2013 ) also find no substantial, massive 
star-forming clouds within the spiral arms. We note that this is also 
seen in the simulations, but only at particular times and locations. 
For example at times 260 Myr and 265 Myr (Figure 3), we see in 
our chosen arm section massive clouds situated in spurs, whilst the 
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Figure 3. The evolution of Cloud380 is shown over a period of 40 Myr. All the panels are shown in the rotating frame of the spiral potential. Clouds are 
indicated by blocks of colour overplotted on the column density plots. The clouds shown are all those which contain at least 20 original cloud particles of 
Clouds 80. Clouds 80 is shown in the middle panel. The fraction of particles which are original cloud particles of CloudS80 are indicated for some of the 
clouds, for clarity fractions are typically only shown for a couple of clouds. The total fraction of the original cloud particles of CloudS80 which are located in 
clouds is also indicated on the top left of each panel. The remainder of this gas lies in the ambient ISM. The locus of the region containing the original cloud 
particles of CloudS80 is indicated by the dashed white lines for each panel (for the middle panel obviously this locus is simply the region in red). 



gas has not yet had time to re-form massive clumps in the spiral 
arm (This again reflects our finding the molecular clouds in spiral 
arms are not 'ridge clouds'). 



relatively small scales (e.g. 10-20 pc) and would require zooming 
in to the plots in Figure 3 to see clearly. 



Stellar feedback contributes to the structure of the gas along 
the spiral arms, and in the inter-arm regions, although clear holes 
relating to multiple supernovae are not easy to see. This is partly 
because we choose a moderate star formation efficiency (for exam- 
ple in pobbs et aL|201l"a| with a star formation efficiency of 20% 
we could clearly see large holes). There is substructure due to feed- 
back within, and around, the clouds in our simulations but it is on 



5.1.3 Cloud formation: evolution at times t < To 

Prior to time To — 20 Myr less than 15 percent of Cloud380 is in 
cloud gas, with the rest in the ambient ISM. In Figure 3 we plot at 



^ Note that we have chosen to plot the cloud components as opaque blocks 
of colour in order to emphasise particular subsets of cloud particles, but that 
this then tends to obscure the detailed structure of the cloud itself. 
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early times (and also late times) the area containing most of the par- 
ticles which find their way into Cloud380. We can see that at time 
To — 20 Myr, most of the gas which eventually forms Cloud380 
lies in the interarm region, and that the interarm gas is generally 
diffuse, rather than lying in clouds. There are a few interarm clouds 
which enter the spiral arm, as seen at 235 Myr, but generally the 
interarm gas is not dense enough or not sufficiently concentrated 
(recall we have a minimum cloud mass of around 10^ Mq) to be- 
come apparent as clouds. 



5.1.4 Does shear or stellar feedback cause Cloud380 to 
disperse ? 

We see that after To = 250 Myr, Cloud380 splits up into a number 
of smaller clouds, and by 270 Myr, a successor of Cloud380 is not 
clearly distinguishable. Since the virial parameter of Cloud380 is 
a = 3, in this instance the cloud is unbound and will naturally ex- 
pand without feedback. However bound regions within Cloud380, 
and other bound > 10^ M© clouds in the simulation (see Figure Al 
in the Appendix and Table 1) cannot be dispersed simply by the 
dynamics of the clouds at To. Also, internal motions which keep 
clouds unbound, likely originate at least in part from feedback and 
travel through the spiral arms (Dobbs et al. 2011a). Overall it is 
likely that Cloud380 splits up overall due to some combination of 
shear and feedback. We argue here that shear is responsible for the 
large scale dispersal of the cloud, but that feedback is also likely to 
break up the densest parts of Cloud380. 

One indication that shear is likely to play a role in the dispersal 
of Cloud380 is simply the morphology of the region at later times, 
and how gas is stretched into a spur. In Figure 4 we show panels 
from a re simulation from 250 Myr where we did not include stellar 
feedback (although in order to do this we also were not able to 
include self gravity). Figure 4 shows panels from the resimulation, 
for the three times of 255, 260 and 265 Myr (equivalent to those in 
Figure 3). The evident similarity of the morphology of the region 
with and without feedback suggests that at least on larger scales, 
it is shear which disrupts the massive clouds and causes them to 
break up into smaller clouds (although because self gravity is not 
included the clouds have dispersed more in Figure 4 compared to 
Figure 3). 

We also consider quantitatively the likely role of shear for 



Cloud380. In Figure 5, we compute the surface density at which 
the tidal force from the galactic potential becomes comparable to 
the self gravity of the gas. More specifically, we compute the criti- 
cal surface density, T^crit such that 

dF 

f^cloud~, — ^critG (2) 

dr 

where dF/dr is the gradient of the force across the cloud due to 
the gravitational potential, per unit mass. For Cloud380 we took 
Tcioud = 85 pc. Figure 5 indicates that at the galactic radius of 
3 kpc, the surface density at which shear becomes important (~ 
90 Mq pc~^) is comparable to the surface density criterion for the 
clouds (100 Mq pc~^). We also show on Figure 5 the timescale 
associated with the shear at different galactic radii. We estimate 
this timescale as the inverse of the Oort's constant A, where 

At R = 3 kpc, this timescale is ^ 30 Myr. As we see in Figure 3, 
this is not too dissimilar for the time for the region to be sheared 
into a spur. Previous authors have attempted to use a Toomre, 
or Toomre-like criterion to determine the relevance of shear in a 
molecular cloud ( Hu nter et al.|1998||Dib et al.|2012| ). However this 
is only applicable for an axisymmetric disc, and there is no equiv- 
alent local criterion. |Dib et al?] (^2012) conclude that shear is not 
important for clouds in the molecular ring, however they used a dif- 
ferent criteria, and obtained slightly lower values of A than found 
in our model. 

Finally we also examined the role of shear and feedback by ex- 
amining the distribution of surface densities within Cloud380 and 
its successors (Figure 5, lower). The most dense regions disappear 
very quickly after 250 Myr, likely due to feedback processes. How- 
ever on longer timescales, of few tens of Myr, the whole distri- 
bution of S is shifted to low surface densities. Again this likely 
corresponds to the whole cloud being sheared out to lower surface 
densities. 

We note that whilst at R ^ 3 kpc, shear appears to be impor- 
tant at disrupting GMCs, at larger radii, this may not be the case. At 
larger radii, Scr^^ is much lower, indicating that dense clouds are 
unlikely to be disrupted by shear. Likewise, the timescale for dis- 
ruptions increases to several tens of Myr. We consider other clouds 
in Section 5.2 
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Figure 5. The top panel shows the critical surface density for shear to be ef- 
fective (calculated according to Eqn 2), plotted against galactic radius (solid 
black line). The shear timescale, calculated as the inverse of Oort's constant 
A is plotted as the blue dashed line. Cloud380 is situated at a radius of ~ 3 
kpc. The lower panel shows the distribution of surface densities in lOpc x 
10 pc boxes over all clouds which contain 20 or more original particles of 
Cloud380, from times of 250 Myr (i.e. for Cloud380) to 270 Myr. There is 
an immediate disappearance of the highest density regions from 250 to 255 
Myr, likely due to feedback. There is then a more gradual decrease across 
all surface densities, likely due to shear. 
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Figure 6. The mass of the continuation of Cloud380 (the largest mass of 
original cloud particles situated in a precursor or successor of Cloud380), 
and several smaller clouds is shown versus time. The time resolution is 
1 Myr between 245 and 255 Myr, and 5 Myr otherwise. 
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5.1.5 Cloud lifetime 

Figure 3 demonstrates that assigning a lifetime to Cloud380 is dif- 
ficult. Cloud380 comprises a number of smaller clouds at earlier 
and later times, and furthermore these clouds do not contain all 
the same particles as Cloud380. So for example at 235 or 270 Myr 
there is no obvious counterpart to Cloud380. We first estimate a 
dynamical lifetime based on the time for the cloud mass to change. 
Figure 6 shows the change in mass of the continuation of Cloud380. 
Recall that the continuation of a cloud is the cloud which contains 
the greatest number of original cloud particles. In Figure 6 we show 
the mass of the original cloud particles in these clouds. Taking our 
definition of the lifetime as the time over which there is an ob- 
ject which has at least half of the same mass of Cloud380 (Section 
4.3.3), then the lifetime of Cloud380 is -20 Myr. 

An alternative lifetime could be based on the time gas remains 
molecular, i.e. based on the chemistry of the cloud rather than the 
flow of gas. In Figure 7 we show the time after 250 Myr for which 
gas exhibits a minimum molecular fraction, for the original cloud 
particles of Cloud380. Because, due to resolution, we insert feed- 



Figure 7. This figure shows the distribution of times the original gas parti- 
cles of Cloud380 remain above a given molecular gas fraction. 



back at comparatively low densities, the most dense gas in our cal- 
culation is only lO's or lOO's cm~^. Gas at the former densities is 
typically not fully molecular (see Dobbs et al. 2008). Gas retains 
a high molecular fraction for only a short time - reflecting that gas 
does not remain very long at densities of >100 cm~^ (Dobbs et al. 
2012), likely due to feedback. Gas retains high molecular fractions 
for up to — 10 Myr. Gas which exhibits molecular fractions > 0.01, 
thus >10 cm""^, retains this fraction for typically — 16 Myr, which 
is more in agreement with the timescales found from dynamical ar- 
guments. The timescale associated with the molecular gas is very 
likely to depend on the surface density of gas in the galaxy and the 
amount of molecular gas. For example, for M51, which is predom- 
inantly molecular in the centre ( Rickard & Palmer|1981| |Scoville| 
|& Young||1983| ) a GMC lifetime based on how long gas remains 
molecular will likely be long, whilst a dynamical lifetime would 
likely be shorter (i.e. arm GMAs/GMCs are not the same clouds as 
inter- arm GMCs seen at some time later). 
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5.1.6 Star formation in Cloud380 

In Figures 8 and 9 we show the star formation which takes in place 
in the gas particles which constitute Cloud380, i.e. the original 
cloud particles of CloudSSO. Although we do not include star par- 
ticles in the simulations, we do record when, where and how many 
stars form each time we insert feedback. Figure 8 (upper panel) 
indicates that star formation is much higher for the gas between 
^240 and ^260 Myr, in agreement with our very approximate 
maximum cloud lifetime of 20 Myr. The total mass of stars formed 
is ^ 6 X 10^ Mq, indicating that the star formation efficiency of 
Clouds 80, measured as the total mass of stars formed (over ^ 20 
Myr) divided by total mass of the gas is around 2.5 %. 

The instantaneous star formation efficiency, although not a 
true efficiency, is an observable measure of the star formation activ- 
ity in a cloud, taken as the mass of young stars present divided by 
the mass of gas. The instantaneous efficiency in nearby clouds (Lu- 
pus, Taurus, Orion, Ophiuchus) is observed to range from < 0.3 % 
to 10 % Puff & Stabler 2006, Merm et al. 2008, M izuno et alT] 
|1995[|Hatchell etal. 2012 ). By comparison, at 250 Myr in our sim- 
ulation, the instantaneous efficiency is ^ 1 %. 

In Figure 9 (top) we show the age distribution of stars formed 
from the original cloud particles of Cloud380 at a time of 250 Myr 
to a maximum age of 20 Myr. Without star particles, we are not 
able to tell when stars decouple from the gas. Instead our results 
represent the maximum ages of stars present in the clouds. Our 
results agree with |Hartmann et al.| ( |2012| ) who find that clouds 
typically have many more very young stars. Our distribution ex- 
tends to higher ages however, as we take a much more massive 
cloud. Tracing the star particles will only make this distribution 
steeper. Changing the feedback scheme (e.g. including a delay) 
could potentially allow more older stars to be present. However 
the timescale for massive star formation is likely to be short, and 
whilst supernovae occur at the end of a massive star's lifetime, stel- 
lar winds are present throughout. 
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Figure 8. The top panel shows the cumulative mass of stars formed versus 
time for Cloud380 (top), which has a mass of 2 x 10^ Mq. The lower 
panel shows the cumulative mass of stars formed versus time for clouds of 
~ 10^ M0. 



5.2 Examples of other clouds 

So far we have concentrated on the evolution of a single cloud. 
Clouds 80. However the behaviour of this cloud may not be typical 
behaviour for clouds in our simulation. In particular, Cloud380 is 
one of the most massive clouds in our simulation. Hence Cloud380 
can only form by the agglomeration of smaller clouds (and ISM 
gas), and can only disperse by splitting into smaller clouds (and 
ISM gas). Thus Cloud380 is effectively at the top of the food chain 
of clouds in our simulation. Smaller clouds can form and disperse 
in a greater number of ways (a simple schematic is shown in Fig- 
ure 10). Smaller clouds can form by the agglomeration of smaller 
clouds, and form instantly as the byproduct of a dispersing larger 
cloud. Similarly smaller clouds can disperse due to feedback and/or 
shear, and disappear instantly by becoming part of a larger cloud. 
Thus the 'cloud lifetime', and evolution of the smaller clouds, can 
be somewhat different from both Cloud380 and usual evolution 
scenarios for molecular clouds (see e.g. [Vazquez- Semadeni et al.| 
|2010j). In this section we consider a number of other examples. 

In Table 1 we show the details of five other clouds we studied 
from the simulation. We tried to select clouds which contained at 
least a few 100 particles, and from a range of environments. We in- 
clude a cloud in the outer part of the disc, where the spiral arms are 
not very strong, and a cloud in the inter-arm region. Unfortunately 
the inter-arm cloud is not so well resolved as inter-arm clouds tend 
to be less massive (see also Columbo et al., in prep.). We do not 



show the detailed evolution of these clouds, but discuss briefly their 
evolution below. The only clouds where the gas does not demon- 
strably form from, or disperse into other smaller clouds, are those 
which are < 10^ M©, which is not surprising as they are towards 
the resolution limit of the simulations, and may well undergo 'arti- 
ficial' formation or dispersal as mentioned in Section 4.3.2. 

5.2.1 Cloud788 

Cloud788 is a 1.7 X 10^ M© cloud in the spiral arm region, which 
over the course of 240-260 Myr, appears to emerge from the dis- 
persal of a more massive (4.3 x 10^ M©) cloud, and then part of 
Cloud788 goes on to form another (2.5 x 10^ M©) cloud. This 
complex evolution is a consequence of Cloud? 8 8 lying in the spi- 
ral arm all this time, where there are merger events, and multiple 
feedback events. The mass of the 'continuation of Cloud788' ver- 
sus time is shown in Figure 6. The lifetime of this cloud is much 
less than Cloud380, being of order a few Myr. 

5.2.2 Cloud877 

Similar to Cloud380, Cloud877 is the highest point in food chain 
during its evolution. However there is only one feedback event dur- 
ing its evolution, suggesting it is shear as it leaves the spiral arm, 
and / or the unbound nature of the cloud, which is causing it to 
disperse. 
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Cloud 


Mass 


No. particles 


a 


Location at 


Nature of cloud 


Nature of cloud 




(105 Mo) 






To = 250 Myr 


evolution 


dispersal 


Cloud380 


20 


6386 


2.9 


Spiral arm (R=3.1 kpc) 


forms from and disperses into smaller clouds 


shear + feedback 


Cloud788 


1.7 


559 


3.7 


Spiral arm (R=4.3 kpc) 


remains of, and progenitor of more massive cloud 


feedback 


Cloud877 


3.1 


999 


1.8 


Spiral arm (R=4.1 kpc) 


forms from and disperses into smaller clouds 


shear + unbound 


Cloud355 


0.96 


305 


3.6 


Inter-arm (R=3.3 kpc) 


remains of more massive CMC 


shear + unbound 


Cloudl59 


2.7 


863 


2.7 


Outer disc (R= 8.3 kpc) 


forms from and disperses into smaller clouds 


unbound 


Cloudll98 


13 


4291 


0.8 


Spiral arm (R=3.4 kpc) 


remains of more massive CMC 


feedback 



Table 1. The different clouds examined in this paper. Outer disc refers to the outer region of the disc (R > 5 kpc) where the spiral pattern is not particularly 
strong. The inner part of the disc is divided into spiral arm and inter-arm. Due to the complexity of determining cloud dispersal, the final column is somewhat 
speculative, but is based on Figure 5, a, the shapes of the clouds, and number / distribution of nearby feedback events (see also text). 
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Figure 9. The age distribution of stars which have formed from the con- 
stituent gas of Cloud380 (top) and Clouds 159 and 355 (lower) are shown. 
The distributions of ages are calculated at 250 Myr. 



5.2.3 Clouds 5 5 

Cloud355 is the remains of a more massive 2.9 x 10^ M©) spi- 
ral arm cloud. Cloud355 is situated in a spur, and simply becomes 
more sheared out with time. There are a couple of feedback events 
between 250 and 260 Myr during the evolution of Cloud355, but 
probably the main mechanism of dispersal is shear. Both Cloud355 
and Cloud877 are unbound, and less massive, so they can more 
quickly disperse. 



5.2.4 Cloudl59 

Cloudl59 appears to form and disperse in a similar manner to 
Cloud380, i.e. from the accumulation of smaller clouds, and disper- 
sal into smaller clouds. However as Cloudl59 lies at a much larger 
radius (8.3 kpc), shear is less likely to be able to disrupt the cloud 
(and would do so on a longer timescale, see Figure 5). In addition, 
there are no feedback events close to when the cloud disperses (un- 
like Cloud380), so feedback is unlikely to be responsible for cloud 
dispersal. Cloudl59 appears relatively elongated (^ 100 pc x20 
pc) and splits into two at the narrowest point of the cloud. In fact, 
Cloudl59 forms from two clouds which adjoin, and then split apart 
again (consequently Cloudl59 only has a very short life, of around 
5 Myr). This suggests the larger scale gas dynamics are determin- 
ing the evolution of Cloud 159. The mass of the continuation of 
Cloudl59 is shown in Figure 5. Similar to the other 10^ M© clouds, 
the lifetime is relatively short. The mass reduces less steeply with 
time, possibly because the cloud is subject to less shear. 



5.2.5 Cloudim 

Cloudl 198 is the remains of a more massive (~ 2 x 10^ M©) cloud. 
Cloudll98 has the longest lifetime, exceeding the duration of the 
40 Myr time period we consider (this is the only such example). 
Unlike Cloud380, even by 270 Myr Cloudl 198 is not particularly 
elongated, and only one relatively small cloud has broken away, 
although the mass of Cloudl 198 is decreasing. This suggests that 
feedback is slowly dispersing the cloud. This cloud is also the most 
bound of those we study individually, and presumably shear and 
large scale motions are not effective mechanisms for dispersing the 
cloud. 



5.2.6 Ages and star formation of clouds 

In Figures 8 and 9 we show the star formation histories for some 
of these lower mass clouds. Figure 8 (lower panel) again shows a 
higher number of stars formed around 250 Myr. There is a peak 
in the star formation rate for between about 5 and 10 Myr, again 
suggesting that a 'lifetime' of these clouds would be several Myr. 
However there are still a number of stars being formed after 255 
Myr, showing that some of the gas is then accreting onto other 
clouds and continuing a small amount of star formation. The mass 
of stars formed from 245 to 255 Myr is - 1500 M© for Cloudl59 
and Cloud788, which gives a relatively low efficiency of < 1 %. 

Figure 9 (lower panel) shows the age distribution of stars in 
Cloudl59 (the distribution is similar for Cloud788 and Cloud877). 
The distribution is more strongly peaked at young ages, compared 
to the more massive Cloud380, and thus more closely resembles 
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Figure 10. A simple schematic of possible cloud evolution is shown. For 
simplicity clouds are represented by solid circles, and regions of ambi- 
ent ISM by dashed ellipses. The evolution of Cloud 4, which forms from 
and disperses into a number of smaller clouds (like Cloud380), is perhaps 
closest to the typical picture of GMC evolution. However the evolution of 
smaller clouds can differ, forming from the dispersal of Cloud 4, or ending 
their lives as a new piece of Cloud 4 (see also Figure l, |Tasker&Tan|2009) . 



[Hartmann et al.| ( [2012| ). Interestingly though, the distribution does 
not go to zero, and there are still a small number of stars with ages 
> 10 Myr. For Cloud355 (Figure 9, lower), which represents the 
remains of a very massive GMC, there is less evidence of a peak 
of very young stars. There is no particular decrease in the fraction 
of stars with age until around 12 Myr, after which there are con- 
siderably fewer stars. We caution again that as we do not include 
any star particles, we cannot say anything about cluster dispersal. 
Nevertheless our results suggest intriguing differences in the stellar 
population of GMCs depending both on their mass, cloud history 
and position in a galaxy. 



5.3 Statistical properties of all > 10^ M© GMCs 

So far we have concentrated on individual examples of clouds. 
Here we extend our analysis to consider collectively all > 10^ M© 
GMCs in our simulation (a total of 61 GMCs). 



40 



30 



6 

£ 20 



10 




10' 











■in DA ^^^^ A 









10^ 
Mass (M ) 



10' 



Figure 11. The lifetime of clouds > 10^ Mq is shown versus their mass. 
Most cloud lifetimes are between 2 and 25 Myr. The clouds with the longest 
lifetimes (> 25 Myr) tend to be the most massive > 10^ M© GMCs. Points 
are shown when applying a critical surface density threshold of 100 (blue 
squares) and 50 (red triangles) Mq pc~^ for the clump finding algorithm. 
For a few of the clouds, a dashed line connects equivalent clouds found with 
the different surface density criteria. The arrows represent lower limits on 
the lifetime of a cloud which exists for longer than the 40 Myr we consider. 



(iv) We define the end of the cloud's life as the last point after 
To when there are at least half of the original cloud particles in a 
successor cloud. 

The time resolution for determining a cloud's lifetime is 1 Myr 
between 245 and 255 Myr, and 5 Myr otherwise. 

In Figure 1 1 , we show our measure of the GMC lifetime ver- 
sus mass (blue points). The GMCs tend to have relatively short life- 
times, the most frequent around 5 Myr, though with a large spread 
from 2 to 40 Myr. The clouds with the longer lifetimes tend to be 
those which are more massive. The lower limit of 2 Myr is a gen- 
uine lower limit on the lifetimes of the clouds (the minimum we 
can resolve is 1 Myr). Clouds can have quite short lifetimes in our 
models because, for example, a cloud that splits into 2 equal parts 
will instantly 'die' in our model. 

We also show in Figure 1 1 points where a lower threshold sur- 
face density (50 M© pc~^) is used to find the clumps. The distribu- 
tion of lifetimes is not dissimilar, though there are relatively more 
10-25 Myr lifetime clouds and fewer clouds with lifetimes ^ 5 
Myr. Equivalent, more massive, clumps found with the lower den- 
sity threshold tend to have longer lifetimes (by a factor of around 2) 
compared to their counterpart clumps found with T^th = 100 Mq 
pc~^, suggesting unsurprisingly that GMAs will likely have longer 
lifetimes than GMCs . 



5. 3. 1 Cloud lifetimes 

We use the same method to calculate the lifetimes as used for 
Cloud380 in Section 5.1.5, but it is stated here more precisely. We 
compute the lifetime of each cloud as follows: 

(i) We follow the continuation of a given cloud by finding the 
constituent precursor or successor cloud at each time frame. 

(ii) We compute the mass of original cloud particles in each pre- 
cursor or successor of our chosen cloud. 

(iii) We define the beginning of the cloud's life as the first point 
before To when at least half of the original cloud particles are situ- 
ated in a precursor cloud. 



5.3.2 Cloud lifetimes and crossing times 

A cloud's lifetime is likely to be related to the velocity dispersion 
of the cloud, and its size, so essentially the crossing time of a cloud. 
|Elmegreen| pOOO| ) suggests that star formation and cloud dispersal 
occurs within 1 or 2 crossing times of a cloud. Of the clouds we 
considered. Cloud 159 undergoes the most straightforward 'termi- 
nation', simply splitting in half at the narrowest point of the cloud 
(see Figure 12). The time for this cloud to disperse is then likely 
to be the length L divided by the relative velocity of the sections 
either side of the part marked L on Figure 12. Taking the velocity 
dispersion of the cloud, a = 5 km s~^, as a proxy for the latter, we 
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Figure 12. The point at which Cloudl59 disperses: at 254 Myr, the suc- 
cessor of Cloudl59 is a single cloud, however at 255 Myr, this cloud has 
split into two. The times are only given to the nearest Myr, but the time 
difference between the two stages is 1 .0 Myr. 
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Figure 13. The actual lifetime of clouds (defined in Section 5.3. 1) is plotted 
against our estimated time given by Equation 4. Outlying clouds are often 
those with particularly elongated, or compact shapes. 



can estimate a timescale oft = L/a. Taking L = 10 pc we obtain 
a time of t =2 Myr. This is a reasonable estimate of the timescale 
to split up the cloud. 

Obviously we cannot study every cloud individually, to deter- 
mine how, and where each cloud splits into multiple pieces. Fur- 
thermore larger clouds such as Cloud380 split up rather by many 
smaller clouds breaking away over a period of time. However we 
can expect that the time for clouds to disperse will depend on their 
dimensions, and their velocity dispersion. We can obtain a simple 
estimate of a cloud's age from the crossing time of a cloud. We 
estimate the lifetime according to 



(4) 



where b is the semi-minor axis of the cloud, and a the velocity 
dispersion of the cloud. Overall, feedback, shear and spiral shocks 
are likely to contribute to a, as discussed earlier, though Equation 4 
makes no attempt to incorporate these explicitly. We include the 



factor 2 because when we select a cloud, it is not at the start of 
its lifetime. We make the simplest assumption that the cloud life is 
symmetric about the point which we observe it. 

In Figure 13 we show the measured lifetime of the GMCs ver- 
sus this estimated lifetime. There is a correlation between the actual 
and estimated lifetimes (the actual lifetime typically within a fac- 
tor of 2 of the estimated lifetime), suggesting that to first order, the 
lifetime is related to the crossing time, however there is quite a lot 
of scatter. Some scatter may simply be because we assumed the 
evolution of the cloud is symmetric about To = 250 Myr, which 
is unlikely to be true for all clouds (see Figure 6). We also inves- 
tigated whether the scatter was related to the aspect ratio, cloud 
porosity, virial parameter or star formation activity and examined 
individual outliers in Figure 13. We found that both the aspect ra- 
tio, and the cloud porosity contributed to the scatter, and outliers 
in Figure 13 were typically either very long, filamentary clouds or 
compact clouds. Thus although we take into account the shape of 
the cloud by using the semi-minor axis in Equation 4, this does not 
satisfactorily represent the complex shape of the cloud, or likewise 
the porosity. Therefore the shape of the cloud appears have a com- 
plex effect on its lifetime, which is difficult to disentangle. We also 
found a weak correlation of lifetimes with galactocentric radius. 
The lifetimes were found to increase slightly with radius, which we 
interpret as being due to the reduced role of shear at larger radii. 

5.3.3 Were (will) most clouds (be) constituent pieces of larger 
GMCs, or are they at the top of the food chain? 

In Section 5.2 we highlighted examples of some ^ 10^ Mo clouds 
that were previously part of larger GMCs. Here we check whether 
GMCs are typically at their most massive when we observe them 
or not. We define 4 types of cloud: 

(i) Clouds which are their most massive at To. 

(ii) Clouds which were part of more massive GMCs prior to To. 

(iii) Clouds which will become part of more massive GMCs af- 
ter To. 

(iv) Clouds which were part of more massive GMCs prior to To, 
and become more massive GMCs again after To. 

We again take our sample of > 10^ M© GMCs at To = 250 Myr 
and look for precursor or successor clouds between To — 20 Myr 
and To + 20 Myr which contain the same particles. For a more 
massive precursor or successor cloud to be classed as significant, 
we denote that it has to contain at least twice as much mass as the 
original cloud. Otherwise, the cloud is assigned to the first cate- 
gory. In Table 2 we show the percentage of different types of cloud. 
Nearly 2/3 of clouds are at least half of their maximum mass at the 
time they are selected (To). 30 % of clouds were previously, or will 
be, part of GMCs which are twice as massive. The remaining few 
are clouds which were previously part of more massive GMCs and 
will again be part of massive GMCs. Thus most GMCs we observe 
are roughly at their maximum mass. However the fraction which 
become, or were previously, part of significantly bigger objects is 
still clearly important. 



6 CONCLUSIONS 

In this paper we have followed the detailed evolution of individ- 
ual GMCs in a global galactic disc simulation. We studied one 
2 X 10^ M© GMC, Cloud380, in particular detail. We showed that 
this GMC was formed from a mixture of clouds (mostly lying in 
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Fraction of clouds at Fraction of clouds previously Fraction of clouds which will form Fraction of clouds which are at 

least 1/2 their maximum mass at Tq part of more massive clouds more massive clouds in the future at a local minimum in mass 

64% 15% 16% 5% 



Table 2. The percentage of clouds which are at their maximum mass, have previously been more massive, or are yet to reach their maximum mass are shown. 
The last column displays the fraction which were previously more massive, and will again be more massive GMCs within 20 Myr. For our purposes, we 
required that the cloud mass had to increase by at least a factor of 2 to avoid being placed in the first column. 



the spiral arms), and ambient ISM (mostly inter arm). This partic- 
ular GMC is effectively at the top of the food chain in terms of the 
evolution of molecular clouds and GMCs. Cloud380 likewise dis- 
perses into smaller clouds (lying in a 'spur' feature) and ambient 
inter arm gas. The tendency of the most massive GMCs to evolve 
into spurs containing GMCs, before more massive clouds start to 
build up in the arm again, can account for recent observations that 
in some sections of spiral arm, the most massive clouds, and the 
most star formation, occurs in spurs ( [Schinnerer et al.,|2013| ). 

More generally, other smaller clouds can have different evolu- 
tionary scenarios. Clouds can be immediately created by the disper- 
sal (or splitting up) of another cloud. Clouds can also immediately 
lose their identity by becoming part of another cloud. Thus the his- 
tory for different clouds will likely be quite varied. 

Because of this complex evolution of GMCs, that the precur- 
sor or successor of a cloud may not be obvious, and that constituent 
material in clouds changes with time, assigning a lifetime to clouds 
is difficult. Here we assign a lifetime as the time over which, if we 
follow the evolution of a given cloud, its mass does not decrease by 
more than half. We further stipulated that this mass had to be the 
same constituent material as when the cloud was selected. We find 
that the minimum cloud lifetime is 2 Myr, but there is a large range 
in cloud lifetimes. Most lifetimes are between ^ 4 and 25 Myr. 
An alternative lifetime could be based on the length of time gas 
remains molecular, though this would mean molecule rich galaxies 
would have long-lived clouds, regardless of whether the constituent 
gas of those clouds remains the same with time. 

We attempted to predict the lifetime of clouds in our simula- 
tion, by computing their crossing time. Although this can provide 
a rough estimate, the complexities of cloud shapes and dynamics 
leads to a lot of scatter, when plotted agains their true lifetimes. 
Clouds which are very filamentary can easily be broken apart, 
whereas those which have small aspect ratios, and little porosity 
cannot be readily split in half. Furthermore as the lifetimes of more 
massive clouds are longer, the velocity field may well change on the 
timescale of the cloud's evolution. We have here looked at shear and 
feedback as possible mechanisms for dispersing clouds. Simply the 
morphology of spurs in galaxies is sufficient to suggest that shear 
is important. We suggest that shear may be responsible at smaller 
galactic radii, and on larger (tens of pc) scales, whereas feedback is 
more important on smaller scales. 

A caveat to these results is that we included instantaneous 
feedback. In some cases, the lifetimes and dynamics of clouds may 
be governed by the larger scale velocity dispersion of the gas, and 
be less dependent on local feedback events. However we leave it to 
future simulations and tests to fully examine the effect of different 
feedback prescriptions. We also plan to study the GMCs in terms 
of CO and HI emission, rather than total density, in a future paper. 

We also examined the star formation history of molecular 
clouds, although we did not include star particles to follow the evo- 
lution of stars. Nevertheless, we found that the age distributions of 
stars are generally heavily skewed towards younger stars, but there 



is some variation according to the history of the cloud (e.g. whether 
it has formed from the dispersal of a more massive cloud). We also 
showed that over the lifetimes estimated for our clouds, the con- 
stituent gas undergoes much more star formation, as expected, and 
the star formation efficiencies of the clouds (as measured by the 
mass of stars formed divided by the mass of the clouds) are of of 
order 1 %. In future work we will include star particles to investi- 
gate star formation in clouds more consistently. 

Finally, we have also examined the properties of GMCs 
formed in our simulation (Appendix). The properties of the GMCs 
(mass function, rotation, virial parameter) are similar to those in 
previous simulations ( |Dobbs et aLl|2011b[ [Tasker & Tan|[2009l 
|Dobbsj2008 ). We do not find much dependence of GMC proper- 
ties on the criteria used for our clump-finding algorithm, or the 
resolution of the simulation. 
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APPENDIX A: PROPERTIES OF GMCS 

In Pigure Al we show the properties of GMCs found in our simu- 
lation. Compared to our previous results, (^Dobbs et al.|2011a|, we 
now extend properties down to 10^ M©. We show results for three 
different implementations of our clump finding algorithm, where 
we take the minimum surface density required to define clumps, 
Et^, as 50, 100 and 200 Mq pc~^. The properties of the clouds 
as shown in Pigure Al appear relatively similar regardless of the 
column density threshold. The slope of the mass spectrum appears 
to become slightly steeper as the threshold increases, the spectrum 
at 12th = 50 M© pc"^ being rather flatter. With 12th = 200 M© 
pc~^ we obtain relatively few clouds, and a less coherent spectrum 
with a few outliers at high masses. There are no low mass clouds 
with Et/i = 200 Mq pc~^ because at least 2 x 10^ M© is required 
in a 10 pc x 10 pc cell (used for our clump-finding algorithm) to 
reach a surface density of 200 M© pc~^. 

The rotation of the clouds is also similar regardless of the 
threshold for the clump-finding algorithm. The main difference is 
simply the cut off in mass, though the values of angular momenta 
also appear slightly lower with T>th — 200 M© pc~^. The val- 
ues of Lz, and relative fractions of retrograde clouds are similar 
to our previous results ( Dobbs 200F,' Dobbs et al.|201 la| ). The frac- 
tion which are retrograde are 39 % (Et^ =50 M© pc"^), 40 % (Et^ 
=100 M© pc"^) and 44 % (T^th =200 M© pc"^). The most massive 
clouds (for any particular Et^) tend not to be retrograde though. 

The virial parameters of the cloud are shown in the third panel 
of Pigure Al. As found in |Dobbs et al.| ( [2011b| ), the majority of 
clouds tend to be unbound. The proportion of bound (a < 1) clouds 
increases from 10% at the lowest threshold to 20% at the highest 
Et/i, which is not surprising as we would expect more dense clouds 
to be bound. 

The number of clouds does obviously differ with the surface 
density threshold chosen. Por the simulations here, this simply re- 
flects the higher fraction of the gas which lies above the lower 
threshold. 
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Figure A2. The column density for the galaxy disc is shown for a lower 
resolution simulation with 1 million particles. 
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Figure Al. The cloud mass function (top), cloud rotation (centre) and virial 
parameters (lower) are shown for clouds selected from our the simulation 
presented in this paper. We tried three different column density thresholds 
to extract the clouds, Et/i=50, 100 and 200 Mq pc~^. In the middle panel, 
the symbols represent Y^th =50 Mq pc~^ (points), =100 Mq pc~^ 
(asterisks) and T^th =200 Mq pc~^ (boxes). For T^th =200 Mq pc~^, 
there are no low mass clouds. This is because we would require at least 
2 X 10^ M0 in a 10 pc X 10 pc cell (used for our clump finding algorithm) 
to reach a surface density of 200 Mq pc~^. 



Al Resolution test 

We also performed a simulation with 1 million particles to see how 
the structure of the disc and properties of clouds change with reso- 
lution. Due to the non-trivial dependence of our feedback prescrip- 
tion on resolution, it is not possible to run an exactly equivalent 
simulation. However, as shown by Figure A2, the structure of the 
disc, and appearance of inter- arm spurs are similar at the different 
resolutions, albeit there is less dense gas, and less obvious dense 
GMCs at the lower resolution. 

We compare the properties of GMCs formed in the 8 and 1 
million particle simulations in Figure A3. In each case we take 
Hth =50 M0 pc~^. We chose this value because taking a higher 
T^th produced a too small number of clouds in the low resolu- 
tion simulation to make meaningful comparisons. We see from Fig- 
ure A3 that there are some differences between the different resolu- 
tions, although we stress again that the simulations will not be ex- 
actly equivalent in the amount of feedback. The mass spectrum for 
the lower resolution simulation (top) is slightly steepe][jThis could 
be related to the finding in the previous section that for stricter cri- 
teria (relative to the resolution and properties of the simulation) the 
spectrum becomes steeper. It is unclear how significant these dif- 
ferences are though. 

The distribution of the specific angular momenta of the clouds 
is very similar in the different resolution simulations. The only dif- 
ference is that there is no obvious truncation in the mass of ret- 
rograde clouds in the lower resolution simulation compared to the 
higher resolution simulation (see also |Dobbs et al.| (^2011a^), it is 
not entirely clear why this is. 

Finally the virial parameters of the clouds are shown in the 
lower parameter of Figure A3. The range of a is similar for 
both resolutions. The clouds appear slightly more bound in the 
lower resolution simulation, probably indicating that the feedback 
is slightly stronger in the higher resolution simulation. The large 
spread in a for the lowest mass clouds appears to be an artefact 
(see also pobbs et al.|2011b| l, though there always appears to be a 



^ The power spectrum also looks more like a curve than a power law. De- 
partures from simple power laws, e.g. multiple slopes or truncated power 
laws are often found to be better fits for GMCs in nearby galaxies ^Hughes] 
|et al.,|2013| . 
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decrease in the spread towards higher masses. The latter may well 
be a selection effect - those clouds which tend to reach the highest 
masses are those which are the most bound. 

This paper has been typeset from a TgX/ I^TgX file prepared by the 
author. 
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Figure A3. The cloud mass function (top), cloud rotation (centre) and virial 
parameters (lower) are shown for clouds from simulations with 1 and 8 
million particles, using J2th=50 Mq pc~^ in our clump-finding algorithm. 
In the middle panel, the results for the 1 million particle simulation are 
represented by asterisks, and those for the 8 million particle simulation with 
dots. 
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